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Electrons incident from a normal metal onto a superconductor are reflected back as holes - a 
process called Andreev reflection m- In a normal metal where the Fermi energy is much larger 
than a typical superconducting gap, the reflected hole retraces the path taken by the incident 
electron. In graphene with ultra low disorder, however, the Fermi energy can be tuned to be smaller 
than the superconducting gap. In this unusual limit, the holes are expected to be reflected specularly 
at the superconductor-graphene interface due to the onset of inter band Andreev processes, where 
the effective mass of the reflected holes change sign SflS]. Here we present measurements of gate 
modulated Andreev reflections across the low disorder van der Waals interface formed between 
graphene and the superconducting NbSe 2 . We find that the conductance across the graphene- 
superconductor interface exhibits a characteristic suppression when the Fermi energy is tuned to 
values smaller than the superconducting gap, a hallmark for the transition between intraband retro- 
and interband specular- Andreev reflections. 

PACS numbers: 78.30.cd,75.47.Lx,75.70.Tj,73.50.Jt 

Andreev reflection (AR) is a process of transferring 
charges from a normal metal (N) into a superconductor 
(SC) W- When an NS interface is biased with an en¬ 
ergy eVns above the Fermi energy ep, an electron can 
only be injected into the SC gap A, if a hole is reflected 
back with an energy of eVns below ep, creating a Cooper 
pair at the Fermi level. For a straight N/SC interface, 
the momentum conservation along the boundary must be 
conserved. Thus the incident angle of an incoming elec¬ 
tron 9inc and the reflected angle of a hole, Oref, have a 
simple relation pe sin Oinc = Ph sin Oref , where Pe and ph 
are the total momentum of the electron and hole, respec¬ 
tively. In the limit, ep ^ A, which holds for a typical 
NS junction, the reflected hole remains in the conduc¬ 
tion band of the normal metal, and therefore necessarily 
carries opposite sign of the mass compared with the elec¬ 
tron. In order to conserve momentum the hole reflects 
back along a path of the incident electron, exhibiting 
nearly perfect retro- AR, with Oref ~ —dine [MS]. 

If however the Fermi energy can be tuned such that 
ep < A, a completely different kind of AR is expected. 

In this case the energy difference, 2eVns, provided by 
the AR process can result in the reflected hole appearing 
in the valence band rather than the conduction band. 

The reflected hole now has the same mass sign as the 
incident electron, and therefore according to momentum 
conservation, travels in the same direction along the in¬ 
terface {Oref > 0 in Fig. I). In the ultimate limit 
of eF=0, the angle of reflection equals to the angle of 
incidence 0ref=0inc - a process called specular reflec¬ 
tion [1H]. The condition for specular AR is satished 


whenever A > eVns > £p (Fig. 1(a) right panel), and 
is therefore predicted to be observable when the normal 
metal in the NS junction consists of a zero-gap semicon¬ 
ductor (ZGS), and the Fermi energy is tuned close to 
the charge neutrality point (CNF) where the conduction 
band and valence band meet. 

Graphene provides an ideal platform to exhibit both 
intra and interband AR. Its two-dimensional nature and 
ZGS properties enable to induce a crossover from inter- 
and intraband AR by tuning ep with the electric field 
effect [MIT]. However, accessing the regime ep ~ A, 
a necessary condition to realize the interband AR, has 
been technically challenging. In typical graphene sam¬ 
ples on a Si 02 substrate, strong potential fluctuations 
up to 5ep ^ 50 meV have been typically observed due 
to the presence of charged impurities [12]. This value is 
much larger than A in typical SC. The recent progress in 
producing suspended [13] and hBN supported graphene 
samples m has now allowed to dramatically reduce these 
fluctuations down to Jej:’ ~ 5 meV [T3]. Despite exper¬ 
imental progress has been steadily made in contacting 
graphene with various superconducting metals |16I420] 
including recent work of the edge contact on graphene 
[nidi], fabricating transparent SC contacts on graphene 
channels with extremely low inhomogeneity has yet to be 
realized. 

In this letter we employ a novel non-invasive approach 
to fabricate N/SC interfaces with an unprecedented en¬ 
ergy resolution close to the CNF. For this purpose we 
electronically couple (see supplementary information (SI) 
for detail method) a high mobility hBN/bilayer graphene 
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FIG. 1. (a) Schematics of the AR processes at a normal metal/SC (left) and a zero-gap semiconductor/SC interface (right). 
An electron from the CB with a total energy of ef + eVns is reflected by a hole with an energy ef — eVns, forming a Cooper 
pair at ef in the SC. For large Ef the hole remains in the CB, resulting in an intraband retro-reflection process. For small 
enough ef the hole undergoes an interband transition into the VB resulting in specular-reflections, (b) Top : color enhanced 
optical image of the hBN/BLG device before the transfer of the NbSe 2 flake (marking its final position). The NbSe 2 and the 
hBN/BLG stack form an electrically coupled overlap region. Bottom left : schematic diagram of the vertical cross-section of the 
vdW stack. Bottom right : sketch of the final device on top of a 300 nm Si02 /Si back gate, (c) Inset : temperature dependence 
of differential resistance dV/dl vs. bias current across the BLG/NbSe 2 interface at I4j=—50 V. Strong variation of dV/dl 
appear below NbSe 2 Tc ~ 7 K. Main : normalized conductance Gi.7k/Giok vs. Vns- The characteristic double-peak line-shape 
arises from ARs in the limit e_F S> A. 


(BLG) device and a 20-100 nm thin NbSe 2 flake. NbSe 2 
is a van der Waals (vdW) SC with a critical temperature 
Tc ^ 7 K and a large A ~ 1.2 meV (Fig. I (b)). 

We use the dry-vdW transfer technique and a current 
annealing method [Ml [26] to achieve ultra-clean, atomi¬ 
cally sharp, and highly transparent Ohmic interfaces be¬ 
tween graphene and NbSe 2 . In this experiment, we chose 
BLG rather than single layer graphene to minimize Ssp 
near the CNP. Owing to the larger density of states of 
BLG near the CNP, smaller Sep can be obtained in BLG 
given the similar level of substrate induced inhomogene¬ 
ity. The devices were fabricated on heavily degenerated 
Si substrates topped with 300 nm Si02, where a back 
gate voltage Vbg is applied to tune ep of the BLG chan¬ 
nel. 

We characterize the electronic transport properties 
across the vdW N/SC junction by measuring differen¬ 
tial resistance dV/dl as a function of the channel cur¬ 
rent Isd (Fig. 1 (b)). Inset of Fig. 1 (c) shows 
typical traces for a high back-gate voltage 14^=—50 V 
that corresponds to a representative condition ep ^ 
A. We observe that the dV/dI{Isd) traces become in¬ 
creasingly non-linear for temperatures T < Tc- In or¬ 
der to resolve these features better, we divide traces 
taken below and above Tc (here at T = 1.7 K and 
T = 10 K) and obtain the normalized differential con¬ 
ductance Gi, 7 /f/Gioic=(dl^/d/ioif)/(dG/d/i. 7 /f). Fig. 1 
(c) shows Gi,tk/Giok as a function of the voltage drop 
across the N/SG junction Vns, estimated by considering 
the BLG channel resistance (SI). The most salient feature 
of the resulting Gi, 7 k/G ioK(Vns) curve is a conductance 
dip around zero bias with two pronounced conductance 


peaks at |I4is| ~ ±1.2 mV. The position of these peaks is 
consistent with the value of A KbSeo, ■ Later, we will show 
that this line shape is characteristic for typical N/SG 
junctions with ep ^ A and can be explained in terms of 
intraband retro-ARs [BHH]. 

The ability to tune ep enables us to investigate the 
characteristic AR signal at the continuous transition from 
large to small ep. Fig. 2 (a) left panel shows the char¬ 
acteristic longitudinal resistance Rxx of the BLG chan¬ 
nel vs. Vbg. An upper bound of dci? < 1 meV was 
estimated from the full-width-half-maximum (FWHM) 
peak of Rxx{Vbg) [SHISl HZ], demonstrating that the 
condition ep < A can be realized in this sample. Fig. 
2(a) right panel displays simultaneous measurements of 
G\, 7 k/G\ok as a function of both the bias current Isd 
and Vbg. For large ep, \Vbg\ > 2 V, Giy^/Gioic(/^d) 
exhibits a similar characteristic behaviour described in 
Fig. 1(c), showing the characteristic conductance dip 
around the zero bias condition Isd = 0. However near 
the CNP, for \Vbg\ < 2 V, the Giyic/Gioic(Ad) traces 
exhibit drastic variations : the conductance dip at zero 
bias turns into a peak and then is suppressed back again 
a,s ep goes through the CNP. 

To explain these experimental findings at the crossover 
from the large to the small ep limit, we develop a theoret¬ 
ical model based on the Bogoliubov-deGennes equations 
and the Blonder-Tinkham-Klapwijk (BTK) formalism for 
the conductance across the BLG/SC interface Gns at 
T=0 K. We also compute the normal conductance Gnn 
for T > Tc 10 K) (see SI for details). Fig. 2 (b) shows 
the resulting normalized conductance Gns/Gnn as a func¬ 
tion of ep/A and eVns/A for comparable energy ranges 
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FIG. 2. (a) Left : R^j; vs. Vbg showing the characteristic 
channel resistance peak at the CNP with a full-width of half- 
maximum of SVbg ~ 0.6 V, suggesting an upper bound of 
5eF <1 meV. Right : color map of normalized differential 
conductance GifkIGwk as a function of bias current Rd 
and the gate voltage Vbg for the inner gap region leV^al < A 
with line-traces of Gi.jk /GioK[lsd) for various fixed Vbg. (b) 
Color map of the theoretically obtained normalized differen¬ 
tial conductance Gns/Gnn vs. ef/^ and Kja/A with line- 
traces of Gns/GnnieVna/A) for various fixed ef/^, showing 
overall agreement with the experimental data. The gate tun- 
ability of BLG allows to continuously probe the conductance 
changes from the large |e_F| S> A to the small |£i?| ~ A limit. 
The zero bias dip line-shape at large ef is continuously trans¬ 
formed into a diagonal cross-like line-shape close to the CNP. 


as the experimental data. The theoretical model demon¬ 
strates qualitatively good agreement with the experimen¬ 
tally obtained normalized conductance in Fig. 2(a). Both 
the broad variation of the normalized conductance in the 
regime of |£i?| A and the rapid variation in the regime 
of IefI ~ A are matched well with our experimental ob¬ 
servation. 

Further quantitative comparison between the experi¬ 
ment and theoretical model for |£f| ~ A can be per¬ 
formed by re-plotting the experimental Gy^kIGi^k map 
as a function of £f and Vns (see SI for the conversion 
scheme) (Fig. 3 (a)). In both graphs, one can iden¬ 
tify four regions of enhanced conductance (colored blue): 
two of them for |£f| > |el4s| and two for |£f| < |eI4s|- 
These regions are separated from each other by connected 
regions of reduced conductance (colored red) that are 
approximately following the dependence |£f| ~ |eI4,j;|, 
forming diagonal lines that are roughly symmetrically 
arranged with respect to ep = eVns — 0 in the con¬ 
ductance maps. Several representative line cuts, show¬ 
ing Gi,-jk/Giqk^sf), clearly exhibit similar features for 
both experimental and theoretical traces (Fig. 3 (b)). 

The various features observed in the conductance map 
can be explained by analyzing the microscopic pro¬ 
cesses for different Sp (Fig. 3 (c)). ARs for SLG and 
BLG involve intervalley processes due to the time rever¬ 
sal symmetry of the backwards motion of the reflected 


hole n mini HE]- Therefore, for ep > eVns, an electron 
in the K-valley of the GB with an energy of ep + eVns 
is reflected as a hole in the AT'-valley of the CB with an 
energy ep — eVns > 0. This intraband AR process gives 
rise to a relatively high conductance, analogous to ARs 
in normal metals in the limits £f S> A. When ep de¬ 
creases, however, the phase space for the reflected hole 
is decreased, resulting in decreasing conductance. This 
effect culminates in a minimum in conductance at the 
condition ep = eR„s, where the hole is at the GNP and 
intraband AR cease to exist. For ep < eVns the hole 
undergoes an interband transition into the VB and the 
conductance again increases. 

This non-monotonic conductance change as a function 
of £f can be quantitatively explained by the existence of 
a critical angle 9^ = sin“^ \/\£f — eIns|/|eF + eI4s| (SI) 
that allows AR process only to happen for electrons that 
are incident at 9inc ^ ^c- In the limit of the intraband 
AR (£f ^ eVns) and the interband AR {ep <R eVns), 
9c ~ nf 2, and thus most electrons undergo AR processes 
resulting in a high conductance. However, for £f ~ eH„s, 
the critical angle approaches zero. If we assume a dis¬ 
tribution of angles for the incident electrons, those out¬ 
side of the critical angle will be reflected normally and 
only those arriving at near normal incidence will undergo 
AR processes, resulting in an overall suppressed conduc¬ 
tance across the junction. For ep=eVns, 6*c=0, and no 
electrons can enter the SG. The resulting conductance 
minima hence mark the exact crossover points between 
intraband and interband ARs. 

We can now connect these processes with the corre¬ 
sponding changes in 9ref (Fig. 3 (d) and (e)). The dis¬ 
persion relation E{pg;,py) of BLG is shown in Fig. 3(d), 
where Px and Py are momenta perpendicular and parallel 
to the N/SC interface, respectively. Gonsidering that the 
parallel momentum py is conserved in the process |4l [5] , 
the excitation spectrum for ARs £ = el4s = \E — ep\ 
can be expressed by a function of only px- Depending 
on the relative size of energy of the incident GB elec¬ 
tron (assuming that it is in the AT-valley) compared to 
£f, we can then relate the reflection process to three dif¬ 
ferent scenarios: (i) Andreev reflected CB holes (normal 
AR), (ii) all CB electrons are normally reflected (no AR), 
and (hi) Andreev reflected VB holes (specular AR), all 
in the AT'-valley. For small eVns ^ ep, Px of the in¬ 
cident electron and the intraband reflected CB hole is 
almost unchanged. Since a CB hole has a negative mass, 
its negative Px and the conserved Py provide an overall 
direction of motion that retraces that of the incident elec¬ 
tron (i.e., (i) above). As Vns grows, px of the reflected 
hole increases, translating into a larger 9ref^ that ulti¬ 
mately becomes 7r/2 when ep = eVns and px = 0. As 
Vns grows further, so that eVns > ep, the incident elec¬ 
tron is reflected as an interband VB hole, which has a 
positive mass and a positive Px that is opposite to the 
sign of Px of the incident electron. While the simulta- 
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FIG. 3. (a) Experimental Gi.jkIGiqk and theoretical GnaIGnn color maps as a function of Vna and ep in the limit |ei?| ~ A. 
A continuous region of lower conductance (red) that is defined for |£f| ~ |eKis| (white dashed lines) subdivides the map into 
four disconnected regions of comparatively high conductances (blue). In the regions where IefI > |eI4is| the ARs are of the 
intraband retro type and in the regions where IefI < leVlisI the ARs are of the interband specular type, (b) Experimental 
Gi.7k/Giqk{sf) and theoretical Gna/Gnn{sF) line-traces demonstrate the evolution of the conductance dips (red arrows) with 
varying Vna- (c) Schematics of the AR process for BLG at the cross-over from intraband to interband ARs. With decreasing 
£f at a fixed Vna the AR hole moves from the GB to the VB. The cross-over point where the hole is reflected onto the GNP is 
defined by eF=eVna- (d) Excitation spectrum e(px)=eVna for a fixed £f < A. With an increasing excitation voltage Vna, the 
momentum of the reflected hole continuously increases from negative to positive values, passing through zero when £F=el4is. 
(e) Schematics of the reflection angles of AR holes in the various energy limits. Starting from perfect intraband retro-reflections 
in the high ef limit, 9 nef continuously increases toward f/ 2 as ef is lowered. At the cross-over point separating intraband 
and interband ARs, EF=eVna, dref exhibits a jump to —7r/2, which eventually results in perfect interband specular reflections 

{6ref = 0inc) wheU £f=0. 


neous sign changes of the mass and of Px have no effect 
on the motion in the z-direction, the sign change of the 
mass abruptly reverses the motion in the y-direction by 
180° due to the conservation of Py, as described in (hi) 
above. For the ultimate limits £f ^ ^ and £f= 0 one 
obtains perfect retro-reflections and specular-reflections 
respectively with varying angles for intermediate energy 
scales and a discontinuous jump at eVns = £f- Summa¬ 
rizing this discussion, the conductance maps in Fig. 3 (a) 
can now be used as phase diagrams separating regions of 
intraband retro-AR, where |£f| > |el4is|, and interband 
specular-AR, where IefI < |eFns|- 

While overall the experimental findings are in qualita¬ 
tive agreement with the theoretical estimates, there are 
certain quantitative mismatches. Notably, the experi¬ 
mental conductance map has an overall higher conduc¬ 


tance with broadened features. In addition, we observe a 
pronounced region of reduced conductance around Vns ~ 
0 for a wide range of Ef that is not present in the the¬ 
oretical map. All these deviations can be attributed to 
imperfections due to a realistically broadened N/SC in¬ 
terface, inelastic scattering at finite temperatures and the 
presence of small potential fluctuations, whose quantita¬ 
tive descriptions go beyond the scope of our simple the¬ 
oretical model. 

In conclusion we have developed unprecedentedly clean 
BLG/NbSe 2 based N/SC junctions that allowed to study 
ARs at low £f- Our observation of gate tunable tran¬ 
sitions between retro intraband and specular interband 
ARs opens a new route for future experiments that could 
employ the gate control of 9 ref, which can be continently 
and independently altered with Vbg and Vns- Most im- 
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portantly our finding help to draw a general picture of 
the exact physical processes underlying ARs. 
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Supplementary Information 


OHMIC ELECTRONIC INTEREACES IN 
HBN/BILAYER GRAPHENE/NBSE 2 
HETEROSTRUCTURES 


To fabricate the hBN/BLG/NbSe 2 hetero-structures 
we pre-fabricate hBN/BLG Hall-bar devices following the 
standard recipe in [14]. We intentionally do not con¬ 
tact one end of the patterned BLG with a Ti/Au elec¬ 
trode, leaving it free for the deposition of NbSe 2 . The 
NbSe 2 flake is then transferred onto the device with the 
vdW transfer technique creating an overlap region with 
the BLG from one side, and four gold electrodes from 
the other side, hence allowing to inject current from 
NbSe 2 directly into BLG. The so created gold/NbSe 2 and 
BLG/NbSe 2 junctions have Ohmic interfacial resistances 
with typical values of 50-500 fl/rm and 1-10 kfl^m re¬ 
spectively (Fig. SI (b)). Fig. SI (a) demonstrates typical 
four-point resistance R vs. temperature T measurements 
of the NbSe 2 flake. 

We estimate the contact resistance Rent between BLG 
and NbSe 2 with the typical four-probe method [521 . By 
injecting a source-drain current /sd=10 nA across the in¬ 
terface we measure the resulting voltage difference Vsd 
across the device (Fig. 1 (b) and inset Fig. SI (c)). 
As one of the voltage probes connects to the pristine 
NbSe 2 and the other to the pristine BLG channel, the 
so measured Vsd signal, in all generality, contains contri¬ 
butions from all three regions - the pristine NbSe 2 , the 
BLG/NbSe 2 interface and the pristine BLG. By neglect¬ 
ing NbSe 2 intrinsic resistivity (as it is in the mSl range 
(Fig. SI (a))) we can estimate the interfacial contact 
resistance through the relation Rent = {Vsd — Vxx)IIsd, 
where we correct for BLGs intrinsic resistivity by mea¬ 
suring Vxx inside the BLG channel separately and sub¬ 
tracting it from Vsd- Overall the so obtained strongly 
back gate dependent Rcnt{Vhg) (Fig. SI (c)) can be also 
directly used for a more accurate estimate of the volt¬ 
age drop across the BLG/NbSe 2 interface Vns{Vbg) = 
Rent ( Vbg ^^sd- 

Here Rcnt{Vbg) has values as low as 400±200 (af¬ 
ter current annealing) at negative Vbg and has asymmetri¬ 
cally higher values at positive Vbg- Overall this behaviour 
is in good agreement with typical metal/graphene junc¬ 
tions [55H5T| and can be explained by pn-junction forma¬ 
tion at the junction edge due to work-function matching 
between the NbSe 2 and the BLG, with the direct conse¬ 
quence that the majority of the current is injected into 
the BLG at the contact edge and not at the contact in¬ 
terface [SO]- Here the work-functions ^NbSe 2 = 5-9 eV 
of NbSe 2 and that of the typically used gold = 4.4 
eV [5S] have comparable values. 


THEORETICAL CONDUCTANCE OF THE 
BLG/SC INTERFACE 

General formulas 

We model the experimental setup by a BLG in the 
normal state in the half plane x > 0 having an induced 
SG gap Ao at X < 0. The proximity to the SC results 
not only in the generation of the gap but also changes 
the chemical potential at x < 0. 

The Bogoliubov-deGennes (BdG) equations describing 
the electron motion in the system can be written as [2] : 


/ n{K^+U (x) -ep A (x) 

^ A*(x) £p - n (^k') - U (x) 

( 1 ) 

where K is the momentum operator and u and v the 
electron and hole eigenfunctions. Assuming that only 
the singlet SC is induced and that only electrons of dif¬ 
ferent valleys form Cooper pairs we write the Hamilto¬ 
nian H (K) of the normal BLG for one valley as a 4 x 4 
matrix [iniMISl] : 



n (K) =hv 


( “ 

ATe-^ 


-t± 


V 0 


ATe^ 

0 

0 

0 


-t± 0 \ 

0 0 

0 ATe-^ 

ATe^ 0 / 


( 2 ) 


where a (K) = aictan (ky/kx)- Here ky is component 
of the momentum K parallel to the interface, kx is the 
x-component of K and v the Fermi velocity. 

Due to the SC and the SC gap A (x) the potential 
U (x) is only finite at x < 0 and these quantities can be 
written in the form : 


f -Do, X < 0 
( 0, X > 0 


(3) 


r Aoe*-^, X < 0 

[ 0, X > 0 


(4) 


The conductance of a SLG/SC junction has first been 
considered by C.W.J. Beenakker mis], while the same 
quantity was calculated for BLG in a subsequent publi¬ 
cation in by T. Ludwig |28j . However, in the latter work 
only the limit of hvtd_ ^ Uq was considered. All this does 
not allow to make a detailed comparison of that theory 
with the present experimental data. 

The sharp change of the potential U (x) and the SC gap 
A (x) at X = 0 is an oversimplification. In the experiment 
the dependence of these quantities on the coordinates 
can be rather smooth. Although this may result in a 
certain overestimation of the transmission amplitude at 
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FIG. SI. (a) Resistance R vs. temperature T of NbSe 2 . The sample undergoes a CDW phase transition below T^dw ~ 33 K 
and becomes SC below NbSe 2 ’s Tc ~ 6.8 K. Both, the “hump”-like feature in R at the CDW transition and the sharpness 
of the SC transition, with a width of ST ~ 0.1 K (inset), demonstrate the high purity of the used NbSe 2 crystal, (b) Top : 
atomic force microscopy image of the atomically flat BLG/NbSe 2 interfacial region forming a sharp N/SC junction. Bottom 
: typical current-biased four-terminal I-V measurements across the BLG/NbSe 2 junction for T > Tc- For all Vj,g the TV 
traces are linear, hence showing Ohmic characteristics, (c) Extracted contact resistance Rent of the BLG/NbSe 2 interface as a 
function of Vbg (with error bars), before (red line) and after (grey) current annealing of the junction with lanneai = 2 mA. The 
current annealing effectively reduces Rent at the same time sharpening the region around the CNP. Inset : schematics of the 
measurement setup. 


the interface, we do not think that the dependence of the 
conductance on ep and eVns can essentially be different 
in a model with a more smooth potential U (x). 

Calculating the conductance one considers scattering 
of particles having the energy e = eVns (counted from ep) 
and moving from right to the left. Here, in accordance 
with the experiment we assume that Uq >> ep, Aq, el4ts, 
with the largest energy in being the coupling energy fivt± 
between the layers hvt± ^ Uq. As the experimental mea¬ 
surements have been presented for e \ Vns\ < Aq, we cal¬ 
culate the differential conductance Gns (e) for jej < Ag. 
In this region, using a formalism developed by C.J.W. 
Beenakker|33j. one can reduce the calculation of the con¬ 
ductance Gns (e) of a N/SC interface to calculation of 
the transmission t (e) and reflection amplitudes r (e) for 
the same interface but in the absence of the SC gap. 

According to this theory the differential conductance 
Gns (£) takes the form : 


Gns (e) = 2GoTr [m (e) (e)] (5) 

where Gq = 2e^ jh and : 

m(e) = ti2 (e) [l - e"^*^r22 (-e)r22 (e)] ^i2i(-^) ( 6 ) 

with : 


(7) 


Eq. (5][7| show that calculating the transmission 
ti 2 (s), t 2 i (e), and reflection ri 2 (e), r 2 i (e) amplitudes 
for right and left moving particles at the interface be¬ 
tween normal metals one obtains the differential conduc¬ 
tance Gns (s) of the N/SC interface. The trace Tr over 
the scattering channels in Eq. <§ reduces in the present 
model to integration over ky. Overall Eqs. ([5j[^ demon¬ 
strate that in order to calculate the differential conduc¬ 
tance Gns, one should simply understand the scattering 
on the interface between two normal metals. 


Experimentally, it is convenient to plot not the con¬ 
ductance Gns {£) itself but the ratio Gns (e) /Gnn (e), 
where Gnn (s) is the conductance of the interface at tem¬ 
peratures exceeding the SC transition temperature Tc. 
Gnn {£,£f) can be easily expressed in terms of the trans¬ 
mission coefficient as : 


(e, £f) = SGoTr |ti 2 (e, e_F)|^ (8) 


At the same time, the conductance G„„(e, £f) given 
by Eq. (|^ is applicable at low temperatures only. In the 
present experiment the temperature at which the normal 
conductance is measured exceeds the energies £p and Ag. 
In this case one should take a thermodynamic average 
{Gnn i£,£F))p using the standard formula : 


P = arccos (£/Ag) 
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{Gnn {S, SF, T))p — 


1 


Gnn (^7 A^) 


d/i 


cosh^ 


(9) 


Eq. Q is used in the present work for the calculation 
of the physical conductance at T > Tc- At the same 
time, the measurements below have been performed 
at sufficiently low temperature, such that one can use the 
formulas derived at T = 0. 

One can see from Eqs. ([^ |^ an important difference 
between the conductance of SLG considered in [3] and 
that of BLG. In the former case the transmission am¬ 
plitude is of order one for any Uq due to the so called 
Klein effect and is equal to a normally incident electron, 
whereas in the latter case the transmission amplitude de- 
cays proportionally to Uq ' at large Uq and can be very 
small. Then, for BLG the conductance Gnn is propor¬ 
tional to \t\^, while Gns, being proportional to |t|^ can 
be much smaller than Gnn- 


In Eq. ( [T^ i ?2 (K) describes the conduction band (up¬ 
per part of the spectrum) and E 4 (K) describes the va¬ 
lence band (lower part of the spectrum). Using the in¬ 
equality fivtj_ ^ Uq we consider only these low lying 
bands. 


Wave functions in the region a; > 0. 

Starting with the Hamiltonian 'H(k), Eq. ([^, and 
using a diagonalization procedure of [35], we can write 
out 4-component vectors u satisfying the equation : 


{H {kx — ep) u = eu (13) 

We consider the case > 0 but the solutions depend 
on the sign ot e + sp and we write them separately for 
e + ep > 0 and e + ep < 0. Eor compact notations we 
introduce the quantity : 


Transmission t (e) and reflection r (e) amplitudes in 
the absence of a SC gap. 

Eigenvalues 


Kq = ^/\£ -I- £p\t_L/hv 
We obtain for the low lying eigenvalues : 


(14) 


Study of the model with the Hamiltonian T-L (K) can 
give results for an arbitrary relation between tp and char¬ 
acteristic energies inside one layer. 

The eigenenergies E (K) of the Hamiltonian T-L (K) in 
Eq. (I) take the form : 


Ei{K) = hv{-tp/2-8{K)), (10) 

E2{Vi) = hv{-tpl2 + 8{K)), 

EQ{Vi) = hv{tp/2 + 8(K)), 

A 4 (K) = hv{tp/2- 8{K)) 

where : 

8{K) = ^{tp/2f +K^ (11) 

Provided that t± exceeds all the other energies, the 
bands of the Hamiltonian ([^ with the spectra Ei (k) 
and Eq (k) are far away from the Fermi surface and their 
contribution into physical quantities can be neglected. 
The eigenenergies of the first two low energy bands take 
in the limit of small kx t_L the following form : 


E 2 {kx) = kv (-tj_/2 -I- 8 {kx)) 


hvK^ 

t± 

hvK^ 


( 12 ) 


£2 — —£p + ~K , 
t± 

It-v 2 

£4 = —£p - K 

t± 


(15) 

(16) 


In order to calculate the wave functions one should 
choose an eigenvalue £ and determine K as a function 
of £. It is clear that constructing plain waves in the 
region of the normal metal Ni and £p > 0 one should 


take the solution (15) for AT at £ > —£p and of (16) 
at e < —£p. In addition, one has a solution for K of 


(16) at e > —£p and of (15) at £ < —£p. However, the 


latter solutions are exponentially growing or decaying as 
functions of x. Nevertheless, they should also be taken 
into account when matching functions at the interface 
because of the deep potential —Uq at a; < 0 where the 
exponential growth can change to plain wave behavior. 

Plain wave solutions at £ -\- £p > 0. In this region we 
have left and right moving electrons from the conduction 
band with the energies from Eq. (15). The solution rtf 
for right moving particles in this region takes the form : 


u 


R 

1 


^ik^x+ikyfj 

COS a 


( KQ/tp \ 

e-“ 

KQjtp 

\ / 


(17) 


while the solution for the left moving particles uf reads 


E 4 {kx) = hv {tp/2 - 8 {kx)) 


tx 
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g —ifca;X+ifc„2/ 

COS a 


/ -Ko/t^ \ 


V 


-Ko/t, 


(18) 




The wave functions uf and U 2 correspond to the eigen¬ 


value E 2 (k) in Eq. (12) and belong to the conduction 


band. They are normalized assuming the current 1 along 
the x-axis for right moving particles and — 1 for left mov¬ 
ing ones. 

Having fixed > 0 we have to express and ky in 
terms of Kq and a. Here we introduce the angle a as : 


Decaying and growing solutions at s + ep < 0. In ad¬ 
dition to the plain waves in Eqs. 18), there are two 
other solutions uf and uf corresponding to the eigen¬ 
value £2 from the conduction band (here as the normal¬ 
ization does not play any role, we omit the pre-factors) 


y< _ ^Ka;X+iKyy 


! -iKojtp \ 

-iK^ltp 

\ / 


and : 


(24) 


kx 'iky - 


(19) 


Then, for £ -|- > 0 we have the following relations 

for the variables a and kx, with a varying in the interval 
—7r/2 < a < 7r/2 : 


kx = Kn cos ( 


ky = Kq sin a 


( 20 ) 


Plain wave solutions at e + ep < 0. For these ener¬ 
gies we have plain waves corresponding to right and left 


moving holes from the valence band in Eq. (16). For 


the right moving holes we obtain the normalized wave 
functions : 


.R 


Jk'^x+ik' y 


j \ 


-Koltp 

V -e“' y 


( 21 ) 


cos a' 

where again the angle en! varies in the interval —p 12 < 
y' < 7r/2. : 


/ iKo/tp \ 




0-1 


V 


iKfiltp 




{ 2 b) 


With the parameters Kx and Py can be written in the 
form : 


Kx = Kq coshy, Ky = Kq sinhy 


(26) 


Decaying and growing solutions at e + ep > 0. In this 
region the growing and decaying wave functions corre¬ 


spond to 64 in Eq. (16) describing the valence band. We 
write the growing itj and decaying functions as : 




/ iKoltp \ 

e^' 

-iKo/tp 

V 


and : 


(27) 


( 22 ) 


Here, the opposite signs in Eq. (201 and Eq. (22) are 


due to the fact that we now consider holes from the va¬ 
lence band instead of electrons from the conduction band. 
In general the current of the right moving holes equals 
to 1, while the current for the left moving holes equals to 
— 1. The solution for the left moving particles takes 
the form : 


^-ik'^x+ik' y 


( -K^/tp \ 

Aol 


2 y/Ko/tp cos a' 


Ko/tp 


V -e 


(23) 


u> = 


/ -iKo/tp \ 
e-^' 
iKo/tp 

V -e^' 


The parameters k'x and kL are : 


k'x = Ko coshy', k' = Ko sinhy' 


(28) 


(29) 


Overall the wave functions written in the region x > 
0 can also be used in the region a: < 0 after shifting 
Ep ^ ep -\-Uo- We denote these functions by adding the 

R L R L 

subscript Uo, thus obtaining ^ 2110 ^ 
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Transmission t 2 i and reflection r 2 i amplitudes. 

Now we calculate the transmission ^21 {s) and reflection 
i'2i (s) amplitudes that match the wave functions written 
on the left and on the right of the interface. Again, we 
should consider the regions e + ep > 0 and e + ep < 0 
separately. 

Region ep > 0, e + ep > 0. The scattering process 
in N2 includes a plane wave incident from the left, 
and another one, reflected from the interface. At 

the same time, region N2 has an additional solution that 
is growing with x (decaying from the interface) with the 
symmetry of uf from Eq. (24). 

After scattering on the interface of Ni one obtains an 
outgoing wave with the structure rtf, Eq. (f?), and a 
decaying wave with the structure , Eq. (28). We de¬ 


scribe the scattering process for e -|- > 0 and match 

these functions at the interface : 


+ ^22 (e) ufc/o + Bu<p^ = t 2 i (e) uf + Cu^, 


(30) 


Eq. (30) presents a system of 4 linear equations. They 


can be simplified using the fact that Uq ^ s,ep and the 
amplitudes ^21 (s) and r 22 (s) can be found denoting : 


Using the same approximation as in (31) and proceed¬ 


ing in the same way as for e + £p one comes to the fol¬ 
lowing results, that again fulfil the Unitarian condition in 


Eq. (35) : 


hi (e) = —2\/ L (e) cos Of' sin a' exp (i$ (a')) (37) 

with the reflection amplitude r 22 (e) : 


^22 (e) = 1 — 2iL (e) cos a' exp (z$ (a')), (38) 


and : 


$ (a') = — arcsin (sin^ a') (39) 

It is well known |34j that the transmission amplitude 
ti 2 (e) is related to the amplitude <21 (e) as : 


ti2 (s) = hi (e) 


(40) 


where no knowledge of the explicit form of the phase 
(5 (e) is necessary for the calculation of the conductances. 


Kn 


|e -l- £p\ 


= L (e) < 1 


we obtain using this approximation : 


(31) 


^21 (e) = 2^L (e) cos a (l -I- sin^ a) exp (—(a)), (32) 
where the angle 4) (a) varies in the interval 0 < 4> (a) < 

p 12. : 


4> (a) = arcsin (sin^ a) . 
We obtain for the reflection coefficient : 


(33) 


Differential conductance of BLG. 


Using Eqs. ([5|^, the calculations of the conductances 
Gns and Gnn are straightforward. From the calculated 
^21 (s) and r 22 (s) one can easily obtain t 2 i (~£) and 
?'22 (—£)■ We concentrate now on the case e > 0, > 0 

where we can distinguish two regions in the parameter 
space, where £ < £p and the reflections are of the retro 
type, and £ > £p where the reflections are of the specular 
type. 


Explicit expressions for the conductance Gns (e) at e < ep 
(retro reflection). 


r 22 (e) = 1 - 2L (e) {l + sin^ a) exp (-■«$ (a)) (34) 

with the Unitarian condition : 


1^21 ie)f + \r22 (^)|^ = 1 


that immediately follows from Eqs. (32 34). 


(35) 


Region ed-e^’ < 0. In the region £ + £p < 0 matching 
the wave functions at a; = 0 results in the equation : 


ufpg + r22 (e) u[p^ + B (e) u<p^ = hi (e) + G (e) u^, 


( 36 ) 


Eqs. (!][§ can be rewritten using the integration over 


ky instead of the trace over the transversal channels. 
However, it is even more convenient to integrate over the 
incident angle a. The angle a corresponds to the energy 
£, while another angle a' corresponds to the energy —e. 
These angles are related to each other by the condition 
that ky is the same in both the cases : 


sin a' Kq (e) 
sin a Kq (—e) 


£p +£ 
£p — £ 


(41) 


resulting in the condition that Andreev reflections are 
only possible for angles |a| < ac, where : 
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a = arcsin 


Sf — s 
Sp ^ 


With : 


(42) 


Qc = arcsin 


e — sp 

S + £p 


(50) 


For $ (a) we use Eq. (43) and obtain for $ (o') : 


$ (a) = arcsin (sin^ a) , 


and : 


(43) 


$ {a') = arcsin 


£ ~\~ Sp \ . 2 

Sin a 

£ — £f 


(51) 


$ {a') = arcsin (sin^ a') = arcsin 


£F -\- £ \ .9 

Sin a 

£F — £ 


(44) 


the conductance Gns can be reduced to the form : 


^ ^ Yg {a, a') cos a 

Gns = ^GqKq (e) / —^^- -yrda 


lo 2\X,(a,a% 


where : 


n (a, a') = 1^21 {e)\^\t2i{-e f = 


(46) 


= 16L (e) L (—e) cos a cos a' (l + sin^ a) (l + sin^ a') 
and: 

Xe (a, a') = ^ [1 - (“£) ^22 (e)] (47) 


Here we can write entering Eq. (46) only up to 


quadratic terms in L (e) : 

|X£(a,a')|^ =sin^/3- (48) 

—2 sin j3 L (e) Vi + sin^ a sin (^d + 4) (a)) + 

+L (—e) \/l + sin^ a' sin (^0 — 4) {a' j) + 

+P {e) (1 + sin^ a) + P {—e) (l + sin^ a') + 

+2L (e) L (—e) \J{^+ sin^ a) (l + sin^ a') x 
X cos (2/3 — 4) {a') + 4) (a)) — 2 sin 4> (a) sin 4> {a') 


Explicit expressions for the conductance Gns (s) at e > ef 
( specular reflection). 

Here the angles a and a' are related to each other as : 


sm a 


Ko{e) 


£ £p 


sina Ko{—£) \ £ — £f 

and the critical angle Oc equals : 


(49) 


Conductance for £ > £f. In this case the conductance 


Gns {£) is determined by Eq. (45). Therefore we obtain : 


YPa,a') = \hi (e)|"|t 2 i (-e)l" = 

= 16L (e) L (—e) cos a (l + sin^ a) cos a' sin^ a' 


(52) 


(45) and as {a, a') is determined by Eq. (47) we get : 


|Xe (a, a')!^ = sin^/3- 


(53) 


—2sin/3 L (e) \/l + sin^ a sin (;0 + 4> (a)) — 

—L (—e) cos a' cos (/3 + 4> (a')) + P (e) (l + sin^ a) + 

+P (—e) cos^ a' + 2L (e) L (—e) cos ol Vl + sin^ a x 
X sin (2/3 + 4> (a) + 4) {a')) — 2 cos (j) (a) sin 4) {a') 


Conductance Gnn of the interface between two normal 
metals. 


The transmission amplitude ti 2 (e) determines the con¬ 
ductance Gnn between two normal metals. At a fixed £f 
the conductance Gnn (e) is given by the following formula 


2 

Gnn {£) = GqKq {£) |ti2 (e) | COS odo (54) 

J — 7r/2 


which with Eq. (32) results in : 


Gnn (e) = GqKq (e) / 4L (e) cos^ o (l -f sin^ a) da 

J —7r/2 

(55) 


= YGoXo(e)L(e) 


The conductance between of the normal metals can 
conveniently normalize the conductance Gns and allows 
to compare with the experimental quantities : 


Gns ( e ) = 


Gns ( e ) 


{Gnn ( e , £F^ T)) j 


(56) 




































7 


It is important to emphasize that, although the not. Therefore, the reduced conductance Gns (e) is not 

conductance Gns (£) is symmetric with respect to the symmetric either, 

inversion e —>■ —£, the normal conductance Gns is 



